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Abstract 

We consider the spreading dynamics of two nested invasion clusters on an infinite tree. 
This model was defined as the chase-escape model by Kordzakhia and it admits a limit 
process, the birth-and-assassination process, previously introduced by Aldous and Krebs. 
On both models, we prove an asymptotic equivalent of the extinction probability near 
criticality. In the subcritical regime, we give a tail bound on the total progeny of the 
preys before extinction. 
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1 Introduction 

The chase-escape process is a stochastic predator-prey dynamics which was studied by Ko- 
rdzakhia [2] on a regular tree. In a previous paper, Aldous and Krebs [1] had introduced 
the birth-and-assassination (BA) process. The latter model can be seen as a natural limit of 
the chase-escape model. In [8] the two models were merged into the rumor scotching process. 
The original motivation of Aldous and Krebs was then to analyze a scaling limit of a queueing 
process with blocking which appeared in database processing, see Tsitsiklis, Papadimitriou 
and Humblet [22]. As pointed in [8], the BA process is also the scaling limit of a rumor 
spreading model which is motivated by network epidemics and dynamic data dissemination 
(see for example, [IT], [5], [18]). 

We may conveniently define the chase-escape processes as a SIR dynamics (see for ex- 
ample p2] or [5] for some background on standard SIR dynamics). This process represents 
the dynamics of a rumor /epidemic spreading on the vertices of a graph along its edges. A 
vertex may be unaware of the rumor /susceptible (S), aware of the rumor and spreading it as 
true/infected (I), or aware of the rumor and trying to scotch it /recovered (R). 

We fix a locally finite connected graph G = (V,E). The chase-escape process is described 
by a Markov process on X = {S, I, R} v . If {u, v} G E, we write u ~ v. For v E V, we also 
define the X — >• X maps I v and R v by : for x = {x u ) u& v, {Iv{x)) u = (Rv(x)) u = x u , if u ^ v 
and (I v (x)) v = I, {Ry{x)) v = R. Let A G (0, 1) be a fixed infection intensity. We then define 
the Markov process with transition rates: 



K(x,I v (x)) = \l(x v = S)J2 1 (*u = I), 
K(x,R v (x)) = l(x v = I)J2 1 ( x u = R), 
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and all other transitions have rate 0. The absorbing states of this process are the states 
without /-vertices or with only / vertices. In this paper, we are interested by the behavior of 
the process when at time there is a non-empty finite set of / and -R-vertices. 

In [13], this model was described as a predator-prey dynamics: each vertex may be empty 
(S), occupied by a prey (I) or occupied by a predator (R). The preys spread on unoccupied 
vertices and predators spread on vertices occupied by preys. If G is the Z rf -lattice and if 
there is no R- vertices, the process is the original Richardson's model [19]. With R- vertices, 
this process is a variant of the two-species Richardson model with prey and predators, see for 
example Haggstrom and Pemantle [12], Kordzakhia and Lalley |15j . 

The chase-escape process differs from the classical SIR dynamics on the transition from 
/ to R: in the classical SIR dynamics, a /-vertex is recovered at rate 1 independently of its 
neighborhood. 

Chase-escape process on a tree If the graph G = T = (V, E) is a rooted tree, the process 
is much simpler to study. We denote by the root of T. For the range of initial conditions of 
interest (non-empty finite set of / and R- vertices), there is no real loss of generality to study 
the chase-escape process on the tree obtained from T by adding a particular vertex, say 
o, connected to the root of the tree. At time 0, vertex o is in state R, the root is in state /, 
while all other vertices are in state S (see figure [1]). We shall denote by X(t) G {S, I, R} v our 
Markov process on the tree T*. Under P^, X is the chase escape process on with infection 
rate A. 




Figure 1: The initial condition : the root is /, o is R, all other vertices are S. 

We say that the Markov process X gets extinct if at some (random) time r < oo, there is 
no /-particle. Otherwise the process is said to survive. We define the probability of extinction 
as 

q T (X) = F X (X gets extinct). 

Obviously, if T is finite then (?t(A) = 1 for any A > 0. Before stating our results, we first need 
to introduce some extra terminology. 

There is a canonical way to represent the vertex set V subset of = U£L N fc 
with N° = and N = {1, 2 • • • }. If k > 1 and v 6 V is at distance k from the root, then 
v = (ii, • • • ,ifc) G V Pi N fc . The genitor of v is • • • ,ik-i)'- it is first vertex on the path 
from v to the root of length k. The offsprings of v are set of vertices who have genitor v. 
They are indexed by ■ ■ ■ , ik, 1), • • • ,(«].,••• , ik,n v ), where n v is the number of offsprings 
of v. The ancestors of v is the set of vertices ■ ■ ■ , i^), < £ < k — 1 with the convention 
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iq = 0. Similarly, the n-th generation offsprings of v are the vertices in Fn N fc + n of the form 
Recall that the upper growth rate d G [1, oo] of a rooted infinite tree T is defined as 



where Vj~ = V n N fe is the set of vertices at distance k from the root and | • | denotes the 
cardinal of a finite set. The lower growth rate is defined similarly with a liminf. When the 
lim inf and the lim sup coincide, this defines the growth rate of the tree. 

For example, for integer d > 1, we define the d-ary tree as the tree where all vertices have 
exactly d offspring^. Obviously, the d-ary tree has growth rate d. More generally, consider 
a realization T of a Galton- Watson tree with mean number of offsprings d G (1, oo). Then, 
the Seneta-Heyde Theorem |21[ [13] implies that, conditioned on T infinite, the growth rate 
of T is a.s. equal to d. For background on random trees and branching processes, we refer to 



For integer n > 1, we define T* n as the rooted tree on V obtained from T by putting an 
edge between all vertices and their n-th generation offsprings. For real d > 1, we say that T 
is a lower d-ary if for any 1 < 6 < d, there exist an integer n > 1 and v G V such that the 
subtree of the descendants of v in T* n contains a [<5 n ]-ary tree. Note that if T is lower d-ary 
then its lower growth rate is at least d. Also, if T is the realization of a Galton- Watson tree 
with mean number of offsprings d G (1, oo) then, conditioned on T infinite, T is a.s. lower 
d-ary (for a proof see Lemma [23] in appendix). 

The first result is an extension of |14l Theorem 1] where it is proved for d-ary trees. It 
describes the phase transition of the event of survival. 

Theorem 1. Let d > 1 and 

Ai = 2d - 1 - 2 v / d(d- 1). 

If < A < Ai and the upper growth rate of T is at most d, then (?t(A) = 1. If \> \\ and T 
is lower d-ary, then < (/r(A) < 1. 

Note that in the classical SIR dynamics, it is easy to check that the critical value of A is 
A = l/(d — 1). Also, for any d > 1, Ai < 1 and, 



The proof of Theorem Q] will follow a strategy parallel to [141 [1] . We employ techniques 
akin to the study the infection process in the Richardson model. They will be based on large 
deviation estimates on the probability that a single vertex is I at time t. 

To our knowledge, there is no known closed form expression for the extinction probability 
<7t(A). O ur nex t result determines an asymptotic equivalent for the probability of survival for 
A close to Ai. Our method does not seem to work on the sole assumption that T has growth 
rate d > 1 and is lower d-ary. We shall assume that T is a realization of a Galton- Watson 
tree with offspring distribution n and 



d = lim sup | Vfc 



l/k 



[SEDJ. 




(i) 




k=l 



It would be more proper to call this tree the complete infinite d-ary tree. 



3 



We consider the annealed probability of extinction: 



g(A) = E'[q T {\)] = P^pf gets extinct) 



where the expectation E'(-) is with respect to the randomness of the tree and P^(-) = E' (Pa(0) 
is the probability measure with respect to the joint randomness of T and X. Note that in the 
specific case d integer and ir(d) = 1, T is the d-aiy tree and the measures ¥' x and Pa coincide. 

Theorem 2. Assume further that the offspring distribution has finite second moment. There 
exist constants cq,ci > such that for all \\ < X < 1, 



Note that the behavior depicted in Theorem [2] contrasts with the classical SIR dynamics, 
where 1 — q(X) is of order (A(d — 1) — 1)+. This result should however be compared to similar 
results in the Brunet-Derrida model of branching random walk killed below a linear barrier, 
see Gantert, Hu and Shi [IT] and also Berard and Gouere [7]. As in this last reference, our 
approach is purely analytic. We will first check that q(X) is related to a second order non-linear 
differential equation. Then, we will rely on comparisons with linear differential equations. A 
similar technique was already used by Brunet and Derrida [9], and notably also in Mueller, 
Mytnik and Quastel [161 section 2]. The parallel with the Brunet-Derrida model of branching 
random walk killed below a linear barrier is that we may think the growing recovery cluster 
as a randomly growing barrier. 

In the case < A < Ai, the process X stops a.s. evolving after some finite r. We define Z 
as the total number of recovered vertices (excluding the vertex o of T^) at time r. It is the 
number of vertices which have been infected before the process reaches its absorbing state. 
We define the annealed parameter: 



The scalar 7(A) can be though as a power-tail exponent of the variable Z under the annealed 
measure P^. In particular, for any < 7 < 7(A), from Markov Inequality, there exists a 
constant c> such that for all t > 1, P' A (Z > t) < ci~ 7 . We define 
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c^e awd-i)) 1 /* < 1 - q(X) < cie 2(d(d-i)) 1/4 



with 



s/X-Xx. 



7(A) = sup {u > : E' X [Z U ] < 00} . 




Theorem 3. (i) For any < A < Ai 




A 2 - 2dX + !-(!- AVA 2 - 2A(2d - lJTT 
2X{d- 1) 




(ii) Let 1 < u < j-n-, A u = u 2 (d - 1) + 2nd + (d — I), and 



A u - y/Al - Au 2 

An — . 

2u 

If X < X u then E f x [Z u ] is finite. If X > X u , E' X [Z U ] is infinite. 
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It is straightforward to check that (i) is equivalent to (ii). This result contrasts with 
classical SIR dynamics. For example, if T is the d-avy tree, for all A < l/(d — 1) there 
exists a constant c > such that E A exp(cZ) < oo. Here, the heavy-tail phenomenon is an 
interesting feature of the chase-escape process. It also appears in the Brunet-Derrida model, 
see Addario-Berry and Broutin pQ, Ai'dekon [2 J and Ai'dekon Hu and Zindy [3j. Note finally 
that 

7(A) ~a;o min ( — — 7T-r,77r ) and 7(A) ~ AtAl 1. 



(d-l)A' 

By recursion, we will also compute the moments of Z. The computation of the first 
moment gives 

Theorem 4. If < A < Ai and A = A 2 - 2X(2d - 1) + 1, then 



(d-l)(l + A + VS) d-1 

Theorem H] implies a surprising right discontinuity of the function A 1— )• E^Z at the critical 
intensity A = Ai: E A Z = 2d/((d — 1)(1 + Ai)) — l/(d — 1) < 00. Again, this discontinuity 
contrasts with what happens in a standard Galton- Watson process near criticality, where for 
< A < l/(d- 1), E' X Z is of order (1 - (d - 1)A) _1 . From Theorem we may fill the gap in 
Theorem Q] in the specific case of a realization of a Galton- Watson tree. 

Corollary 5. Let T be a Galton-Watson tree with mean number of offsprings d. Then a.s. 
<?t(Ai) = 1. 

The method of proofs of Theorems [30 will be parallel to arguments in [8] on the birth- 
and-assassination process. 



The birth-and-assassination process We now turn to the BA process. It is a scaling 
limit in d — >■ 00 of the chase-escape process on the d-eay tree when A is rescaled in X/d. 

Informally, the process can be described as follows. We start from a root vertex that 
produces offsprings according to a Poisson process of rate A. Each offspring in turn produces 
children according to independent Poisson processes and so on. The children of the root are 
said to belong to the first generation and their children to the second generation and so forth. 
Independently, the root vertex is at risk at time and dies after a random time D that 
is exponentially distributed with mean. Its offsprings become at risk after time D and the 
process continues in the next generations. In other words, conditioned on D , each offspring 
has a independent exponentially distributed life time with mean one. We now make precise 
the above description. 

As above, W = U^i N fc denotes the set of finite k— tuples of positive integers (with 
iV° = 0). Elements from this set are used to index the offspring in the BA process. Let 
S N-*, be a family of independent Poisson processes with common arrival rate A; 
these will be used to define the offsprings. Let {D v },v £ N*, be a family of independent, 
identically distributed (iid) exponential random variables with mean 1; we use them to assign 
the lifetime for the appropriate offspring. The families and {D v } are independent. The 
process starts at time with only the root, indexed by 0. This produces offspring at the 
arrival times determined by H that enter the system with indices (I), (2), • • • according to 
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their birth order. Each new vertex v, immediately begins producing offspring determined by 
the arrival times of H„. The offspring of v are indexed (v, 1), (v, 2), • • • also according to birth 
order. The root is at risk at time 0. It continues to produce offspring until time T = D , 
when it dies. Let k > and let v = (m, • • • , re/c-i, n^), v' = (m, rtjt-i) denote a vertex and 
its genitor. When a particle v' dies (at time T v i), the particle v then becomes at risk; it in 
turn continues to produce offspring until time T v = T v i + D v , when it dies (see figure [2]) . 

The BA process can be equivalently described as a Markov process X(t) on {S, I, R} Nf , 
where a particle/vertex in state S is not yet born, a particle in state / is alive and a particle 
in state R is dead. A particle is at risk if it is in state I and its genitor is in state R. We use 
the same notation than above : under W\, the process X{t) has infection rate A > 0, q{\) is 
the probability of extinction and so on. 




* 



Figure 2: Illustration of the birth-and-assassination process, living particles are in red, dead 
particles in blue, particles at risk are encircled. 

The following result from |4] describes the phase transition on the probability of survival 
as a function of A. 

Theorem 6. Consider the BA process with rate A > 0. //AG [0, 1/4], then q{\) = 1, while 
tfX>\,0<q(\) < 1. 

The critical case A = 1/4 was established in [8]. Note also that the threshold A = 1/4 is 
consistent with ([1]). 

Our final result is the analog of Theorem [2j 

Theorem 7. Consider the BA process and assume that A > 1/4. There exist constants 
co,ci > such that for all 1/4 < A < 1, 

cowV^ 1 < 1 - q(\) < cio^V^ -1 , 

with 
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Note that the analog of Theorems [3]|3] was already performed in [8] . The remainder of 
the paper is organized as follows. In section [2j we prove Theorem [TJ In section [31 we prove 
Theorem [2j In section [H we prove Theorems OH Finally, in section [5j we extend our previous 
arguments to the BA process and prove Theorem The proof of Theorem [7] can be read 
independently. It is slightly simpler than the proof of Theorem [2j 



2 Proof of Theorem [T] 

We define the set recovered and infected vertices as R(t) = {v G V : X v (t) = R} and I(t) = 
{v G V : X v (t) = I}. The set R(t) being non-decreasing, we may define R(oo) = Ut>oR(t) 
and Z = \R(oo)\ G NU{oo}. Note that also a.s. R(oo) = {v G V : 3t > 0,X v (t) = I}, in 
words, R(oo) is the set of vertices which have been infected at some time. 



2.1 Subcritical regime 

We fix < A < Ai. In order to prove that <7t(A) = 1> it is sufficient to prove that E^Z < oo. 
Let Vk be the set of vertices of V which are at distance k from the root 0. The time for 
infection of v G V n is the sum of n i.i.d. Exp(A) variables say £i, • • • , £ n . The time for recovery 
of the ancestor of v is the sum of n i.i.d. Exp(l) variables, say D\, • • • , D n , independent of 
(£i> " " " j Cn)- Let vq, ■ ■ ■ ,v n be the ancestor line of v. vq = and v n = v. The vertex v will 
have been infected if for all 1 < m < n, v m -\ has infected v m before being recovered. We 
thus find 

(m m \ / n n \ 

i=l i=l / Vi=l i=l / 

The Chernov bound gives for any < 9 < A, 

(n n \ ( / n n \l 

i=l j=l / I \i=l i=l / J 

A \ n / 1 xn 



.A-ey v i + i 

where we have used the independence of all variables at the last line. Now, the above expres- 
sion is minimized for 6 = (A — l)/2. We find 

4A 



(v G R(oo)) < 

Also, by assumption, 



(A + 1)2 

V n \ < (d + o(l)) n . 



It follows that 



»ev n=0 
It is now straightforward to check that 



< 1, 



(A + 1)2 

if A < Ai . This concludes the first part of the proof. 
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2.2 Supercritical regime 



We now fix A > Ai. We should prove that <7t(A) < 1. The monotony of the process in A 
implies that we may assume without generality that Ai < A < 1. For 5 > 0, we define the 
function g$ by, for all x > 0, 

g s {x) = --log(i) +--log(-) -2-log(<5) 

X \X J X \x J 

1 + A , [x 2 \ 

+ log — - 2. 



x \X8 

Taking derivative, the minimum of gs is reached at c = (1 + A)/2. We deduce easily the 
following property of the function <jy. 

Lemma 8. If \% < A < 1, min x> Q gd(x) < 0. 

By Lemma [8l using continuity, we deduce that there exist c > and 1 < 6 < d such that 

g s (c) < 0. 

In the remainder of the proof, we fix such pair (c, 5) . 



Construction of a nested branching process. We fix an integer m > 1 that we will be 
completely specified later on. We assume that m is large enough such that T* m contains a 
[~<5 m ]-ary subtree. We denote by T' this subtree and by x G V its root. For integer k > 0, we 
define as the set of vertices of generation k in T'. Note that by assumption 

Mi = rn fc - 

We may assume that the generation of x in T is larger than m. We denote a(x) G V the m-th 
ancestor of x in T. For z G and fe > 1, we denote by a(z) G V^._ 1 its ancestor in T". For 
example, if z G V/, a(z) = x. 

The chase-escape process can be constructed thanks to i.i.d. Exp(A) variables (£v)veV an d 
independent i.i.d. Exp(l) variables {D v ) v& y. The variable £ v (resp. D v ) is the time by which 
v £ V will be infected (resp. recovered) once its ancestor is infected (resp. recovered). 

We now start a branching process as follows. We set x to be the root of the process, 
5o = {x}. We define recursively the offsprings of the fe-th generation as the set Sk of vertices 
z G satisfying the following three conditions : 

1. the vertex a(z) G V^_ 1 belongs to Sk-i] 

2- X^gtt & — ~ where 7r is the set of the vertices on the path from a{z) to z (excluding 

a(z)); 

^- X^ie7r-^« — ? where ir is the set of the vertices on the path from a(a(z)) to a(z) 
(excluding a(z)). 

Thus for z G V£, such that its ancestor a(z) G we have that 

(m \ / m \ 
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Notice that by construction, the number of offsprings of z ^ z' in S^-i are identically 
distributed and they are independent if a{z) ^ a(z'). It follows that the process forms a multi- 
type branching process with a finite number of types (one for each possible configuration of 
Si). In the next paragraph, we will check that this branching process is supercritical, i.e. we 
will prove that 

M = Y, Pa(* e S x ) > l. 

v<=V{ 

It implies that with positive probability, the branching process does not die out (see Athreya 
and Ney [6j chapter 5]). 

We first notice the following. Assume that at some time t > 0, the vertex x becomes 
infected and that a{x) is still infected. Note that the existence of such finite time t > has 
positive probability. Let us denote by E such event. We set to = t and, for integer k > 1, 

m 

tk — tk-i H • 

c 

By construction, if E holds and z £ then, at time t^, z and a(z) are both infected. It 
follows that in order to conclude the proof of Theorem [T] it is sufficient to prove that M > 1. 



The nested branching process is supercritical. We need a standard large deviation 
estimate. We define 

J(x) = x — logx — 1. 

The next lemma is an immediate consequence of Cramer's Theorem for exponential variables 
(see [10, §2.2.1]). 

Lemma 9. Let (£j)i>i, be i.i.d. Exp(A) variables. For any a > 1/A, we have that 



while, for any a < 1/A, 



liminf — logP ( > am) > -J(Xa), 

m->-oo m \ ^— ' / 

\i=l / 

liminf — log P f < am] > -J(Xa). 



Note that the bounds of Lemma [U] hold for all a > (even if they are sharp only for the 
above ranges). We may now estimate the terms in (|2|). We have from Lemma [9] that 



\i=l 

and 



^ & < ^ j > exp |— mJ + o(m)| 



(^2 Di - ~ ^ - ex p|- mJ +o(m)|. 
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Thus we obtain a lower bound on the mean number of offspring in the first generation to be 

AT = £> A (*€Si) 

zev{ 

> \5 m ]exp{-m( j(-) + j(-) +o{m) 



c J \c 
> exp (-mg$(c) + o(m)) , 

where gg (■) is as defined in ([2]). If m was chosen large enough, we have that M > 1 and hence 
that the branching process is supercritical. Therefore with positive probability, this branching 
process does not die out. This proves the theorem. ■ 



3 Proof of Theorem [2] 

3.1 Differential equation for the survival probability 

We first determine a differential equation associated to the probability of extinction. Under 
P' x , define Q\(t) to be the extinction probability given that the root is recovered at time 
t > so that 

poo 

?(A) = / QxWe-'dt (2) 
Jo 

and Q x (fi) = 1- 

Now, in T, the offsprings of the root are {1, • • • , N}, where N has distribution tt. The 
root infects each of its offspring after an independent exponential variable with intensity A. 
Let {(,i}i<i<N be the infection times. Note that in T, the subtrees generated by each of the 
offsprings of the root are iid copies of T. Hence, if for integer i with 1 < £j < D , we define 
X 1 as the subprocess on vertices (iW) n V with ancestors i. Conditioned on D = t, on N 
and (£i)i<j<jv, the processes (X 1 ) are independent chase-escape processes conditioned on the 
root becomes at risk at time t — £j (where we say that a /-vertex is at risk if its genitor is in 
state R). 

For the process X to get extinct, all the processes X % must get extinct. So finally, we get 

\{ (1(& >t) + < t)Q x (t - & + A)) 

l<i<N 

where (Di),i > 1, are independent exponential variables with parameter 1. Consider the 
generating function of tt 

oo 
k=0 

Recall that ip is strictly increasing and convex on [0, 1] and = EiV = d. We find, for any 



Qx(t) = K 
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t > 0, 



Q x (t) = i)(e- xt + \J^ e~ Xx J Q x {t-x + s)e- s dsdx^ 
+ Ae- A *y e Xx J Q x (x + s)e- s dsdx) 



Tp(e- Xt + Xe~ xt J e {x+1)x J™ Q x (s)e- S dsdx^ . 



Performing the change of variable 

x(t)=if>- 1 (Q x (t))€[Q,l], (3) 

yields to 

pt f'QO 

x(t) = e~ xt + Xe~ xt / e (A+1)x / ip{x{s))e- s dsdx. (4) 

Jo Jx 

We multiply the above expression by e xt and differentiate once, it gives 

roc 

e xt (Xx(t) + x'(t)) = Ae( A+1 )* / ^(x(s))e- s ds, (5) 



Now, multiplying the above expression by e and differentiating once again, we find that 

x(t) satisfies the differential equation 

x" - (1 - A)£c' + ip(x) = (6) 

with 

<p(x) = Xip(x) — Xx. 

This non-linear ordinary differential equation is not easy to solve. However, in the neigh- 
borhood of A = Ai it is possible to obtain an asymptotic expansion as explained below. 

3.2 A fixed point equation 

We define p E [0, 1) as the extinction probability in the Galton- Watson tree: 

P = i/>(p)- 

We note that ip is convex, <p(l) = <p(p) = 0, (p is negative on (p, 1) and it is increasing in a 
neighborhood of 1, p'(l) = X(d — 1) > 0. 

Let T~L be the set of non- increasing functions / : M+ —> R+ such that /(0) = 1, lim^oo f(t) = 
p. The next lemma is an easy consequence of the monotony of the process. 

Lemma 10. For any X > X±, the function x(-) defined by ([3]) is in %. 

Proof: The monotony of the process implies that x(t) is non-increasing. We may thus define 
a = linn.-j.oo x(t). Using the continuity of ip yields to 

ft roo ft 

Xe~~ xt / e (A+1) * / ijj{x(s))e~ s dsdx = Xe~ xt / e A ^(a)(l + o(l))dx, 
Jo Jx Jo 
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This last integral being divergent as t — > oo, we deduce that 



Xe~ xt / e (A+1)x / tf;(x{s))e- s dsdx = ^(a) + o(l). 



From (j4j, we get that a = ^>(a) which implies that a G {p, 1}. Note however that the inequality 
a < 1 follows from Theorem Q] which implies that q(\) < 1. ■ 



From now on in this section, we fix a small u > and we assume that 

Ai < A < 1 - u. (7) 
We define the map A : U ->• R + ) defined by 

A(y)(t) = e~ xt + Xe~ xt e^ +1 ^ x / V(2/(s))e" s dsdx. (8) 

J ix 

Since 

IHIl°°([o,i],r + ) = max |V(x)| = 1 
xe[o,i] 

it is indeed straightforward to check that A(y) is in A{y){t) is bounded by 1. 

Note also that y = 1 is a solution of the fixed point equation 

y = A(y). 

By ([5]), we find that the function x defined by ([3]) satisfies also the fixed point x = A{x). In 
the sequel, we are going to analyze the non trivial fixed points of A. 

Let x G T~L such that x = Ax. Then x ^ 1. By induction, it follows easily that t t— >■ 
is twice differentiable. In particular, from the above argument, x satisfies ([6]) and we are 
looking for a specific non- negative solution of ([6]) with x(0) = 1. To characterize completely 
this solution, it would be enough to compute x'(0) (which is necessary negative since x(0) = 1, 
x'(0) = corresponds to the trivial solution x = 1). We will perform this in the next subsection 
in an asymptotic regime. We first gives some basic properties which follow from the phase 
diagram of the ODE ©. 

Notice that if at some time t > 0, x'(t) = then from ([6]) and p < x(t) < 1, we deduce 
that x"(t) > 0. In particular, x'(s) > for all s G (t,t + 5) for some 5 > 0. This contradicts 
that x(-) is non-increasing. Also, from ([5|) and A < 1, we find that for all t > 0, 

- 1 < x'{t) < 0. (9) 

We define X(t) = (x(t),x'(t)) and 

F(xi,x 2 ) = (x 2 , (1 - A)x 2 - (f(xi)) 

so that 

X' = F(X). (10) 

We define the trajectory <3? = {X(t) : t > 0}. Recall that p = lim^oo x{t). Also, since for all 
t > 0, X(i)i = F(X(t)) 1 < 0, $ is the graph of a differentiable function / : (p, 1] (-1,0) 
with /(l) = x'(0) < 0, 

$ = {(*,/(*)) :ae(M]}. 
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Figure 3: Illustration of the phase portrait. In blue, the curve T, in red the curve 



Moreover 

We notice that on the curve 

r = {(x u x 2 ) G [p, 1] x [-1,0] : (1 - A)x 2 = (f(xx)} 

the second coordinate of F vanishes (see figure [3j). The next lemma shows that our function 
(x(t),x'(t)) cannot cross T near its origin (x(0), x'(0)). 

Lemma 11. There exists 5 > such that for any A as in ( [7)1 the following holds. Let x G T~L 
such that x = Ax. If (x(t), x' (t)) G T for some t > 0, then \x'(t)\ > 5 and x(t) < 1 — 5. 

Proof: Consider the larger s such that (s,/(s)) G V. Recall that f(l) = x'(0) < 0. Thus, on 
(s, 1], (s, f(s)) is below T and it follows that / is increasing. 

We define a G (p, 1) as the point where the function <p(x) reaches its minimum. Since 
(p is decreasing on [p,a], it follows that s G [a, 1]. We set f(s) = —5 and assume that 
< 5 < \(p(a)\/(l — A), so that s G (a, 1). We will prove by contradiction that 5 cannot 
be arbitrarily small. This will conclude the proof of the lemma. Indeed, by construction 
(1 — A)<5 = —<p(s) and (p has a continuous inverse in a neighborhood of 1 : as 5 goes to 0, 
s = s(5) goes to 1. 

To this end, we fix any (3 G (a, 1) and set b = —(p(/3) > 0. Assume that 5 is small enough 
so that f3 < s. Consider the solution Y(t) = (y(t),y' (t)) of the ODE (|1Q[) with initial condition 
y(0) = (p, -5). The trajectory of Y(t) is denoted by 4> = : t > 0}. We define the set 

r+ = {(x 1 ,x 2 ) G [a, 1) x [-5, 0) : (1 - X)x 2 > (p{x x )}. On T+, < and F(x) 2 > 0. It 

follows that the trajectories $ and <!> exit T + either on its left side {(a, x 2 ),x 2 G [—5, 0]} or its 
upper side {(xi,0),xi G [a, 1)}. From Q, we deduce that $ exit T + on the left side. Since 
<I> and <1> cannot intersect and <i> is on the left side of $ in r+, we deduce that necessarily, $ 
also exits T + on the left side. We now check that if 5 is too small, this is not possible. 

Define r > as the exit time of Y(t) from r+. If < t < r, using that ip is increasing 
on [a,/3], we find y"(t) > —(1 — X)5 + b > 6/2 for 5 small enough. We deduce that for all 
t G [0, t), y'{t) > (b/2)t - 5 and y(t) > {b/4)t 2 -St + p. Set t e = 25 /b. We deduce that r < t e 
and, if 5 is small enough, y(r) > (6/4)tg — 5t e + /3 = —5 2 /b + /3 > a. Hence, if 6 is too small, 
must exit T + on the upper side. This is not possible. ■ 
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3.3 Comparison of second order differential equations 



It is possible to compare the trajectories of solutions of second order ODE by using the phase 
diagram. For two functions <pi,ip2 ° n [0, 1], we write ip\ < <p 2 if for all t G [0, 1], <p\(t) < ipzfy). 

Lemma 12. Let 5 > be as in Lemma \11[ Let x G H such that x = Ax. Let dp be a 

Lip <shitz- continuous function and y be solution of y" — (1 — X)y' + <p(y) = with y(0) = 1, 
y'(0) < 0. We define the exit times 

T = M{t > : (y(t),y'(t)) £ (0, 1] x (-1, 0)}, 

T_ = w£{t > : y'{t) < -1} and T + = inf{t > : (1 - X)y'(t) = (p(y(t)), y{t) > 1 - 6}. 

(i) IfT+ = T < oo and(p>Lp then y'(0) > x'(0). 

(ii) IfT- = T < oo and <p < <p then x'(0) > y'(0). 

Proof: Let us start with the hypothesis of (P. The proof is by contradiction : we assume 
that y'(0) < x'(0). We set Y(t) = {y(t),y'(t)) and G( yi ,y 2 ) = (y 2 , (1 - A)y 2 - <p(yi)). Define 
the trajectories $ = {X(t) : t > 0}, and for < r < T, $(t) = {Y(t) : < t < r}. Recall 
that $ is the graph of an increasing function / given by (lllj) . 

Similarly, if y G (0, 1] x (—1,0), G(y)\ = y 2 < 0. Thus there exits a differentiable function 
5 : (y(T), 1] — ► (—1,0) such that 

$(T) = {(s,g(s)): S e(y(T),l}}, 

with 

,'( S ) = 1-A ^ 



9{s) 

Now, the assumption y'(0) < x'(0) < reads g(l) < f(l) < 0. 

By assumption, there exists 1 — 5 < s < 1 such that (s,g(s)) G {(^1,^2) : (1 — X)x 2 = 
(p(xi),x 2 < 1 — 5}. Hence, by Lemma [11] and the intermediate value Theorem, there exists a 
largest s < s\ < 1 such that the curves intersect : g(s\) = /(si) and g(s) < f(s) on (si, 1]. 

Note however, that from (jlip . if s G (si, 1], 

9 '( S ) = 1-A-M >1 _,_^£l = m 

/(a) 

Hence, integrating over (si,l] we find g(l) — g(s\) = g(l) — f(s\) > /(l) — /(si). This 
contradicts 5(1) < /(l). We have proved ([i|). The proof of ((n]) is identical and is omitted. ■ 



3.4 Proof of Theorem M 

We first linearize ([B]) in the neighborhood of Ai. 
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Step one : linearization from below. We have ip(l) = 0, <p'(l) = \(d — 1) > 0, and from 
the convexity of <p, 

tp(a) > X(d- 1)(«- 1). (12) 
We take Ai < A < 1 and consider the linearized differential equation 

y"-(l-AV + A(d-l)(y-l) = 0. (13) 

The solutions of this differential equation are 

(l-A)t (l-A)t 

y(t) = 1 + a sin(u;i)e 2 +6cos(a;t)e 2 ; 

where 

w = -\/-A 2 + 2(2d-l)A-l = c(A) v / A- Ai, 

and 

c(A) = ^(2d-l) + 2 v / d(d-l)-A = (d(d - 1)) 1/4 + 0(|A - Ail). 
We use this ODE to bound from below a^(O). 
Lemma 13. For all Ai < A < 1, 

x' A (0) > -ei^e-^^l + <3(u; 2 )). 
Proof: Let a < 0, 6 = (1 — A)/2, and consider the function 

y (t) = 1 + osin(u;t)e fe *. 

We have y(0) = 1, y'(0) = auj, 

y'(t) = ae bt (uj cos(wt) + 6sm(wt)), 

and 

y"(t) = ae w (26wcos(wt) + (6 2 - u 2 ) sin(wt)). 

Define 

m 7T 1 / 26w \ 7T 2 ^. 2n 

i = arctan -5 = = + G>(w ). 

On the interval (0, T), < and y"(T) = 0. Thus the function y'(t) is decreasing on [0, T] 
and 

y '(r) = e-fae^(w + 0(w 3 )). 
Hence, we may choose a such that y'(T) = — 1 with 

a = -u~ x el e~f(l + 0(u 2 )). 
It remains to use (1121) with Lemma P 
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Step two : linearization from above. For < r/ < min(l, c(A)/(d — 1)), we define 

£ = (l _ V )X + f?Ai < A, 
and the Lipschitz-continuous function 

<p(s) = max (<p(s),t(d - l)(s - 1)) . 

In particular 

<p<(p. (14) 

We define the linear differential equation 

y" - (1 - \)y> + £(d - l)(y - 1) = 0. (15) 

The solutions of (fT5|) are 

y(t) = 1 + asin(o; t)e 2 + bcos(u)t)e 2 5 

with 

J = i^-A2 + 2(2d - 1)A - 1 - Md ~ 1)(A - Ai) = uJl- V ^~^\ 

A careful choice of a, will lead to the following upper bound. In the sequel, o(l) denotes a 
function which goes to as oj goes to 0. 

Lemma 14. // ir has finite second moment, then there exists a constant cq > such that for 
all Ai < A < 1, 

x' x (0) < -co^e'^^il + o(l)). 

Proof: We set 



b = and k = * 1 — — - - . 

2 V C W 

We parametrize in terms of k, so that 

£ = A — (1 — k 2 )u> 2 and a;' = kuj. (16) 
For a < 0, we look at the solution 

y(t) = 1 + asm(ujKt)e bt . 

We have y(0) = 1, y'(0) = okuj. 

y'(t) = ae bt {ujK cos (unt) + bsin(ujKt)). 
We repeat the argument of Lemma [T3l On the interval [0, T], y"(t) > and y"(T) = 0, where 

7T 1 , / A 7T 2 2 n 



T = arctan — — = + 0(u 

ujk ojk V ir — oj K I ujk b 
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and the O(-) is uniform over all n > 1/2. The function y'(t) is increasing on [0,T] and 

y '(T) = e~ 2 ae^ {lok + 0(lo 3 )). 

Now, we have tp(s) < £{d — l)(s — 1) for all s £ [1 - a, 1] with 

-£{d - l)a = ip(l - a) = X{ip{l -a)-l + a). 

If 7r has finite second moment then, from Abel's Theorem, ip" is continuous on [0, 1]. Also 
from Jensen's inequality, ip"(l) > d{d — 1) > 0. We expand tp in a neighborhood of 1, as 
u — > 0, it yields to, 

* - ^(A-4(i + o(i)) 



<(1)A 



(i-^K(i + o(i)), 



where o(l) is uniform over all < k < 1. In particular, for all u small enough, a < 5 with 5 
as in Lemma [TTJ Also, from f)15[) . for t = T, since y"{T) = 0, we have 

y(T) - 1 1 - A 26 



y'(T) £(d-l) 
We may choose a such that y(T) = 1 — <r by setting 

~£(d—l)e~^,^ ,.. 2 ^(d-l) 2 (1 - k 2 )w 
a = -ae 2 1 + 0{u 2 )) = -e 2 ±— —- ! - e ^ '— 1 + o 1 . 

26 Xlp {1)0 K 

By construction, with this choice of a, we have (1 — A)y'(T) = cp(y(T)). Now, in the domain 
1 — c < y < 1 the non-linear differential equation y" — (1 — A)y' + (p{y) obviously coincides 
with (fTBI) . Thus, using (fT4l) and Lemma CEJ it yields to 

*'(0) < y'{0) = -e 2 ^|e-^(l - k 2 )u 2 {1 + (1)). 

Taking finally k = 1 — uj/{nb) and using (fT6l) give the statement. ■ 



Step three : end of proof. We may now complete the proof of Theorem [2j We start 
with the left hand side inequality. We first note that, by Lemma [TTJ x'{t) is decreasing on 
the interval [0, to] where to is the time where (x(to), a/(to)) £ T = {{x\,X2) G [p, 1] x [—1,0] : 
(1 — A)x2 = ip{xi)}. Moreover by Lemma [TTT we find x(to) < 1 — 5. However, by ([9]), we have 

x(t) > 1 - t. 

Hence to > S. Then, by construction, on the interval [0,to], 

x{t) < 1 + x'(0)t = 1 - \x'(0)\t. 
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Since ip{x) < x on [p, 1], it follows from ([2]) that the survival probability may be lower 
bounded as 

poo poo 

l-g(A)= / (1 - V(x(t)))e~*(it > / (1 - x{t))e- l dt 
Jo Jo 

> / |x'(0)|te~*dt 
Jo 



rS 

> \x'(0)\ / te^dt. 
Jo 



It remains to use Lemma [TH and we obtain the left hand side of Theorem [2j 

We turn to the right hand side inequality. For X = (x\,X2) £ [p, 1] x (— oo,0), define 
G{X) = (x2, (1 — A)x2). From the definition of F in (jlOp . we have, component- wise, for any 
X G [p, 1} x (-oo,0), 

F(X) > G{X). 

Note also that G is monotone : if component- wise X > Y then G(X) > G(l"). It follows 
easily that if X(0) = Y(0), X' = F(X) and Y' = G(Y) then component-wise 

X(t) > Y(t). 

Looking at the solution of y" — (1 — X)y' = such that y(0) = 1 and y'(0) = x'(0), we get that 

x(t) yi + x'ioXe^-^-i). 

We deduce from ([2])-([3|) and the convexity of ip that, 

poo 

q(\) = / ^{xit^e^dt 
Jo 

roo 

> / V(l + z'(0)(e (1 - A)t -l)) e -*dt 

> y°° (l + ^OXe^'-r/Je-*^ 

> l + ^'(0)/A. 

We finally apply Lemma [13] and this concludes the proof of Theorem [2j ■ 



4 Proofs of Theorems 
4.1 Proof of Theorem [4] 

For the tree T^, without counting o, let be the total number of recovered vertices when 
the process reach its absorbing state given that the root dies at time t. By definition, if D is 
an independent exponential variable with mean 1, then 

Z = Y(D), 
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where the symbol = stands for distributional equality. 

In T, we denote the offsprings of the root by {1, • • • , N}. The random variable N has 
distribution it. The root infects each of its offspring after an independent exponential variable 
with intensity A. Note that in T, the subtrees generated by each of the offsprings of the root 
are iid copies of T. Hence, the recursive structure of the tree T leads to the following equality 
in distribution 

TV 

Y(t) i l + <t)Y i (t-& + A)- (17) 

where (£i)ieN are iid exponential variables with intensity A, (y)i<i</v and (A)i<i</v are 
independent copies of Y and D respectively. Note that since all variables are non-negative, 
there is no issue with the case Y{t) = +oo in the above recursive distributional equation 
(RDE). The RDE ()17p is the cornerstone of the argument. 
We start with a lemma 

Lemma 15. Let t > and u>\, ifE' x [Z u ] < oo then E' x [Y(t) u ] < oo. 

Proof: Since Z = Y(D), from Fubini's Theorem, E' X [Z U ] = f™ E' A [Y '(t) u ]e - * dt . Therefore 
E,' x [Y(t) u ] < oo for almost all t > 0. Note however that since t H> Y(t) is monotone for the 
stochastic domination, it implies that E^[Y(£) U ] < oo for all t > 0. ■ 



Now, assume that E' A Z < oo. We may then take expectation in (fT7|) : 

ft roo 

E' A y(t) = l + d / / E A [y(t-a; + s)]e-^sAe~ Ax, dx. 
io io 

Let fi(t) = WyY(t), it satisfies the integral equation, for all t > 0, 

/*£ /'OO 

f 1 (t) = l + \de- xt e {x+l)x h{s)e' s dsdx. (18) 
Jo ii 

Multiplying by e A * and taking derivative, we get: 

poo 

{flit) + Xh(t))e xt = Xe xt + \de^ / /i( S )e-^ S . 



Then, multiplying by e~^ A+1 ^, taking the derivative a second time and then re-multiplying 
by e , we obtain: /{'(£) — (1 — X)f[(t) — A/i = —A — Xdfi(t). So, finally, /i solves a linear 
ordinary differential equation of the second order 

x" — (1 — X)x' + X(d — l)x = —A, 

with initial condition /i(0) = 1. We get that 

h{t) = x{t) 1 



d-V 

where x(t) solves the ordinary differential equation 

x" -(l-X)x' + X(d-l)x = 0. (19) 
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with x(0) = d/(d — 1). The discriminant of the polynomial X 2 — (1 — X)X + \(d — 1) = is 

A = A 2 -2A(2d-l) + l. 
If < A < Ai, the discriminant is positive. The roots of the polynomial are 

a= 1 ~ X -^ and P = ' - " + ^. (20) 
2 P 2 v ; 

The solutions of (fl~9j) are 
for some constant a. 

Similarly, if A = Ai, then a = /3 = l — d+ d(d — 1) and the solutions of f)19|) are 

For < A < Ai, we check easily that the functions x(-) with a > are the nonnegative 
solutions of the integral equation (|18p . 

It remains to prove that if < A < Ai then EZ < oo and 

Indeed, we would get EZ = / °° EF(i)e -t <ii = -^ij^ - ^zj as stated in Theorem H To this 
end, we consider T n as being a Galton- Watson tree with offspring distribution ir stopped at 
generation n. We denote by Y^ n \t) the total number of recovered particles given that the 
root dies at time t. We have Y^(t) = 1 and for all n > 0, as in RDE (fTTl) . 

N 

y(n+l) (t) d 1 + Y^l^< t )Y} n \t-^ + D l ), 
i=l 

where 1$ , and are independent copies of y( n ) and D respectively. Since EN < oo, for 
all n > 0, f[ n \t) = EY^ n \t) < oo. Taking expectation, we have for all t > 0, 

/f +1) (f) = 1 + Ade- Ai f e^ x+1 > r f[ n \s)e- s dsdx = Hf[ n) )(t), (23) 

JO Jx 

where <3? is the mapping 

ft poo 

<S>:g^l + Xde- Xt / e {x+l)x / g(s)e- s dsdx. 

Jo Jx 

It is easy to check that <1> is indeed a mapping from H\ to Hi, where Hi is the set of non- 
decreasing functions g : [0,oo) — > [l,oo) such that sup t>0 g(t)e~ at < oo. Now, from what 
precedes the function 

, , , de at - 1 

m = — ■ 
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is a fixed point of <E>. 

We notice that fi(t) < h(t). Due to the monotony of the mapping $ we get /-[ = 
&(fi°^) < = ^- By recursion, it follows f[ n ^ < h. The monotone convergence Theorem 

In) 

implies that flit) = limn-^oo f\ (t) exists and is bounded by h{t). Therefore /i solves the 
integral equation (fTSj) and is equal to x— l/{d— 1) where x is given by ([2T]) (or (f22|) if A = Ai) 
for some a > 0. However, from what precedes, we get x(t) < h(t) + l/(d — 1) and the only 
possibility is a = and fi(t) = h(t). 

This concludes the proof of Theorem [H ■ 



4.2 Proof of Theorem [3] for integer moments 

For < A < Ai, we define 

_ = A 2 -2rfA + l-(l-A)y/A = 

7 2A(d-l) a' 1 ' 

The key property of 7(A) is that (1 — X)ua — \{d — 1) — n 2 a 2 > if and only if 1 < u < 7(A). 
We also note that if u > 1, u < 7 is equivalent to A £ (0, A M ). We first state an important 
lemma. Let 1 < u < 7, we define H u , the set of measurable functions h : [0, 00) — > [0, 00) 
such that sup 4>0 h(t)e~ uat < 00. Let C > 0, we define the mapping from % u to 

* : h ^ Ce uat + A(ie- A * / e {x+1)x \ h{s)e~ s dsdx. 

JO Jx 

In order to check that ^> is indeed a mapping from T~L U to % n , we use the fact that if 1 < u < 
7 = P/a. then ua < (3 < 1. 

Lemma 16. Lei 1 < it < 7 and / 6 % u snc/i i/tai / < ^f(f)- Then for all t > 0, 

(ua + X)(l-ua)e uat 



f(t) < c 



(1 — A)na — X(d — 1) — n 2 a 



Proof: We set go = f and for > 1, we define = ^(^-i). First, since 1 < u < 7 then 
(ua + A)(l — ua) > Xd. We use the formula for all u > such that na < 1: 



t roo \(„aut „—\t\ 



Xe- Xt I e^ x+1 > I e aus e- s dsdx = / (e °" /" \ . (25) 
' ' (na + A)(l - na) 



Jx 

uat 



We deduce easily that if g (t) < Le then 

9l (t) = *(g )(t) < Ce uat + ^ -(e uat - e~ xt ) < L,e uat , 

[ua + A)(l — na) 

with Li = C + ( Ma+ ^) A ('f_ MQ! ) ■ By recursion, we obtain that limsup fc (^(t) < L 00 e uat , with 
Loo = C(ua + A)(l - na)/((l - A)na - A(d - 1) - n 2 a 2 ) < 00. 

We may now conclude the proof. Notice that \I/ is monotone: if for all t > 0, /ii(i) > /i2(i) 
then for all t > 0, \P(/ii)(t) > ^(/^(i). Hence, by recursion, from the assumption / < ^(/) = 
<7i, we deduce that for all integer k > 1, / < a*;. It remains to take the limit in k. ■ 
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Now, let p be an integer, and define f p (t) = E^[Y(t) p ]. The main result of this subsection 
is the following lemma. 

Lemma 17. Let 1 < p < 7^, if X £ (0, X p ), then f p is finite and there exists a constant C p 
such that for all t > 

/„(*) < C p e pat . 

Proof: In §4.1} we have computed f p for p = 1 and found fi(t) = {de at — l)/(d — 1). Let 
p > 2 and assume now that the statement of Lemma [T7] holds for q = 1, • • • ,p — 1. Let n > 0, 
y( K )(i) = min(Y(t), k) and let < s t denote the stochastic domination, from RDE (|17p . we have 

iV 

Y^{t) < st 1 + < t)Y{ K \t - & + A). (26) 

i=l 

Recall the multinomial formula 

n \ v 

71 1 -P 1 . . . „Pn 



5> = E L.",K-«e 

\i=l / Pl,-,Pn VPl 

where the summation is taken over ro-tuples of integers that sum up to p. Taking power p in 
the above stochastic inequality and expanding brutally, we thus get 

N+1 ) II (>=° + < t)Y} K \t - & + Ar) , 

where the summation is taken over N + 1-tuples of integers that sum up to p. Now we define 

= E' A [Y^(tf] = E' A [min(y(t), k)*] . 
Taking expectation and using independence leads to 

°° — / 4- 1 \ n 

fi?\t) < 5>(n) U Il(^=o+^>i E A n<t)Y^(t-t + DT 

n=0 Pi,-,Pn+i VPl---Pn+l/ , =1 V 

* e*<»> e C;.;:,J 

n=0 pi,-, Pn+1 ^ rn+i; 



n / rt pod \ 

x II (*Pi=o + \ t >i^- Xt j o J f^(s)e~ s dsdxj . (27) 

Consider a n + 1-tuple that sums up to p such that for all i = 1, ■ ■ ■ ,n, pi < p — 1, 
YH=iPi = Q < P and Y17=i\i>^ = 771 — P- From the recursive hypothesis and ([25]) . with 

L 1 = maxi< fc <p_i C k , we get 



n / rt poo \ 

n ^=o + l Pl >iAe- A * / / /W( S )e- s ^x < [] 

i=l V J ° Jx ' i: Pi > 



C Pi \e a ^ 



(A + ma)(l - PiOt) 

"•Pi — 

< L™e aqt e~^= lln ( 1 ~ Pia \ 
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Now recall that | ln(l — y) + y\ < 2 (i- y ) ^ or V e -0- Since Y17=iPi — 9 ' we 

II^=o + l Pl >iAe- A *y e^ x+1 > J f Pi (s)e- s dsdxj < L^e^e aq+ ^=^) 



< L 2 e apt . 

Then, grouping together all such n + 1-tuples, from (|27p we deduce 

00 / ft POO \ 

4 K) (i) < Y j ir(n)((n + l) p L 2 e apt + n\e- xt e {x+1)x f^(s)e- s dsdx) 
n=o V Jo Jx J 

f't f'OO 

< L 2 e apt + Xde~ xt / e {x+l)x f^ K) {s)e~ s dsdx, (28) 



where we have used the hypothesis p < j n . We may then apply Lemma [T6j there exists a 
constant C p such that 

f[ K \t)<C p e apt . 

The monotone convergence Theorem implies that fp(t) — lim n _ fl (t) exists and is bounded 
by C p e apt . The recursion is complete. ■ 

4.3 Proof of Theorem [3] : lower bound on 7(A) 
To prove Theorem [3l we shall prove two statements 

If E' X [Z U ] < 00 then u < 7, (29) 

If 1< u < min(7, 7?r ) then E' X [Z U ] < 00. (30) 

In this paragraph, we prove (|30p . The argument is a refinement of the argument in §4.21 
Let K > and let fi K \t) = E / A [min(y(t), k) u ], we have the following lemma. 

Lemma 18. Ifl < u < min(7, there exists a constant C u > such that for allt>0 and 

K > 0, 

f't /*CO 

ft ] (t)<C u e uat + \de- xt e {x+l)x f^\s)e- s dsdx. 







Proof: The lemma is already proved if u is an integer in (|28p . The general case extends of the 
same argument. We write u = p + v with v £ (0, 1) and integer p > 1. We use the inequality, 
for all 2/j > 0, 1 < i < n, 

(t*)'*t E ( " k + " n tf. (3D 

\i=l / i=l pi,-,p n VFi Fn/ l<j<nj& 

where the summation is taken over n-tuples of integers that sum up to p (which follows from 
the inequality (Y2Ui) v — YlVi an d the multinomial formula). Then from (|26p . we get the 
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stochastic domination 



Y^(tr < st £ £ f^t 1 ) (v=o+wi(ei<^ (K) (i-^+Ar +u ) 

x J] (l P;= o + l P; >i%<^ W (t-^ + Ar), (32) 

l<j<N,j^i 

where the summation is taken over N + 1-tuples of integers that sum up to p. From Lemma 
[T71 there exists C such that for all 1 < q < p, f q {t) < C q e qat . Note also, by Jensen inequality, 

that for all 1 < q < p — 1, f q+v (t) < f p (t) p < C p e^ q+V ' . The same argument (with p 
replaced by u) which led to (|28p in the proof of Lemma [T71 leads to the result. ■ 

Statement (|30p is a consequence of Lemma [16] and Lemma [THJ Indeed, by Lemma [TBI 
for all t > 0, /l K) (t) < Cie uat for some positive constant C u independent of k. From the 
monotone convergence Theorem, we deduce that, for all t > 0, f u {t) < C u e uat . However from 

Z = Y(D), we find 



E' X 7 U 



f'OO poo 

K Z U = / /„(i)e-*dt< / C u e UQ *e"*^. 
Then, statement (f30l) follows from no < 1. 



4.4 Proof of Theorem [3] : upper bound on 7(A) 

In this paragraph, we prove statement (|29p . This will conclude the proof of Theorem [3l Let 
u > 1, we assume that E' A [Z U ] < 00 we need to show that A < X u . Without loss of generality 
we can assume that A < Ai. From Lemma [15] and (|17p . we get 

/TV \ u 

f u (t) = E' x [Y(t) u ] = E' A f 1 + £l(6 < t)y,(t - & + A) J • 

Taking expectation and using the inequality (x + y) M > x u + y", for all positive x and y, we 
get: 



> l + \de- xt e ^ +1 ^ x f u (s) e - s dsdx. 

JO Jx 



(33) 



From Jensen's Inequality, f u (t) > fi(t) u > e uat . Note that the integral e aus e~ s ds is finite 
if and only if u < a . Suppose now that 7 < u < a . We use the fact: if u > 7 then 
u 2 a 2 — (1 — \)ua + X(d — 1) > 0, to deduce that there exist < e < A such that 

u 2 a 2 - (1 - X)ua + X(d - 1) > ed. (34) 

We define A = A - e, A(e) = (1 - A) 2 - 4(Ad - A). Note that A(0) = A. Since A < Ai, for e 
small enough, A is non-negative, we may then consider the real roots of X 2 — (1 — X)X + \d— A: 

1-A- V 7 ! ~ i-A+\/A 
a(e) = and /3(e) = . 
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Again, for e = 0, 5(0) = a and /3(0) = f3. Hence, since u > 7 = /3/a, by continuity, for e 
small enough, 

ua > /3. (35) 
We compute a lower bound from ([33]) as follows: 

ft roo rt /'OO 

fu(t) > 1 + ede- xt / e( A+1 > / f u (s)e~ s dsdx + Ade" A * / e ( A+1 > / f u (s) e - s dsdx 

JO Jx JO Jx 

ft /'OO /*£ /'OO 

> 1 + ede- xt / e ( A+1 J x / e uas e~ s dsdx + Ade~ A * / e ( A+1 ^ / f u (s) e - s dsdx 

Jo Jx Jo Jx 

ft fOQ 

> 1 + L(e uat - e- xt ) + \de~ xt / e^ x+l > / f u (s) e - s dsdx, 

Jo Jx 

with 

£ = 7^ rr, -r > 0- 

(A + ua)[l — ua) 

We consider the mapping f : /i i-> 1+L(e uat — e~ xt )+\ e x f u (s)e~ s dsdx. \P is monotone: 
if for all t > 0, hi(t) > h 2 {t) then for all t > 0, > #(fr 2 )(*)- Since, for all t > 0, 

fu(t) > ^(f u )(t) > 1, we deduce by iteration that there exists a function /i such that h = 
^(h) > 1. As in ^4.11 solving /i = ^(/i) is simple, taking twice the derivative, we get, 



uat 



h" - (1 - \)ti + (Ad - A)/j = -A - L(A + ua)(l - ua)e 

Therefore, h = ae at + be^ 1 — e(u 2 a 2 — (1 — X)ua + \{d — l))~ l e uat for some constant a and b. 
From (j35|) the leading term as t goes to infinity is equal to — e(u 2 a 2 — (1 — \)ua+\(d— l))~ 1 e uat . 
However from (134p . — e(u 2 a 2 — (1 — X)ua + A(d — 1)) _1 < and it contradicts the assumption 
that h(t) > 1 for all t > 0. Therefore we have proved that u < 7. ■ 



5 Proof of Theorem [7] 

5.1 Differential equation for the survival probability 

The proof of Theorem [7] is essentially identical to the proof of Theorem [2j We first determine 
the differential equation associated to the probability of extinction. Define Q\(t) to be the 
extinction probability given that the root dies at time t > so that 



/(A) 



Q A (t)e~*«ft 



(36) 



and Q\(0) = 1. Let {£i}i>i be the arrival times of H with < £1 < £2 < • • • • For integer i 
with 1 < £j < D , we define 2% as the subprocess on particles iW with ancestors i For the 
process B to get extinct, all the processes Bi must get extinct. Conditioned on H , and on the 
root to die at time D = t, the evolutions of the (Bi) then become independent. Moreover, 
on this conditioning, Bi is a birth-and-assassination process conditioned on their root to be 
at risk at time t — £j. Hence, we get 



Qx(t)=E > 



Y[Q\(t-& + Di 



Y[Qx(€i + Di) 



,<t 
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where {£i}i>i is a Poisson point process of intensity A and (-Dj), i > 1, independent exponential 
variables with parameter 1. Using Levy-Khinchin formula, we deduce 

Qx(t) = exp J (EQ A (x + D x ) - l)dxj 

= exp fx J J (Q\(x + s) - l)e~ s dsdx^j . 

So finally, for any t > 0, 

Q\(t) = exp f-Xt + A J e x J Q x (s)e- S dsdx^ . (37) 
Performing the change of variable 

x{t) = -In Q x (t) (38) 

and differentiating ()37[) once yields 

x'(t) = e* / V9(x(s))e _s (is, (39) 



where y(y) = A(l — e y ). Now, multiplying the above expression by e * and differentiating 
once again, we find that x(t) satisfies the differential equation 

x" - x' + <p(x) = 0. (40) 



5.2 A fixed point equation 

Let T> be the set of increasing Lipschitz-continuous functions ip on M + such that <p(0) = 0. 
T>o is the subset of functions in T> such that ||</?||oo < 00 • Let T~L be the set of measurable 
functions / : M+ — >• such that /(0) =0 and for any o > 0, 

lim e~ as f{s) = 0. 

s— >-oo 

For ip £ Vq, we define the map T v : H — )• % defined by 

T^(y)(t) = e x ip(y(s))e- s dsdx. (41) 

Since ||^||oo < oo, it is indeed straightforward to check that T v (y) is indeed an element of 
T~L (Ttp(y)(t) it is bounded by H^H^i). Note also that y = is a solution of the fixed point 
equation 

V = T v {y). 

If ip(s) = A(l — e _s ), then using (|37l) we find that the function x defined by (j38j) satisfies also 
the fixed point x = T v (x). In the sequel, we are going to analyze the non trivial fixed points 
of T v with ip G D. 

Let x such that and x ^ 0. By induction, it follows easily that t i-> x(t) is 

twice differentiable. In particular, since x(s) > 0, x'{t) > and the function x : R+ — )■ R+ 
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is non-decreasing. Moreover, by assumption there exists to > such that x(to) > 0. By 
monotony, we deduce that x(t) > for all t > to. Then, using again (|39p . we find that for all 
t > 0, 

< x'(t) < Halloo. (42) 

From the above argument, we are looking for a specific non-negative solution of (|40p which 
satisfies x(0) = 0. To characterizes completely this solution, it would be enough to compute 
x'(0) (which is necessary positive since x(0) = x'(0) = corresponds to the trivial solution 
x = 0). We first gives some basic properties. We define X(t) = (x(t), x' (t)) so that 

X' = F(X) (43) 

with F((xi, x 2 )) = (x 2 ,x 2 — ip(xi)). We also introduce the set 

A = {(x 1 ,x 2 ) G : tp(xi) < X 2 < IMloo}- 

Lemma 19. Let if G Vq and x £ H such that x = T v x and x'(0) > 0. Then x satisfies ([4*0|) 
and for all t > 0, 

X{t) € A. 

Moreover 

lim x'(t) = Halloo- 

t— >oo 

Proof: Define the trajectory $ = {X(t) G R\ : t > 0}. Since for all t > 0, Xit)^ = 
F(X{t)) x > 0, $ is the graph of a differentiable function / : [0, S) -)• R + with /(0) = x'(0) > 0: 

$ = {(sJ(s)):se[0,S)}, 

with 5 = Hindoo x(t) £ (0, oo]. Moreover 

fl(g) _ F((s,f(s))) 2 _ f(8)-<p(8) 

The graph of the function ip is the curve L = {(s, (f(s)) : s G and the set 

A' = {(x 1 ,x 2 ) G [0,oo) 2 : x 2 < ip{xx)} 

is the set of points below L. Assume that the first statement of the lemma does not hold. 
Then by (|42p and the intermediate value Theorem, the curves L and <3? intersect. Then the 
exists so > such that 

/(s ) = (f(s ). 

From flE}, /'(s ) = while ip'(s ) > 0. It follows that (s,/(s)) G A' for all s G (s ,si) for 
some si > so- By induction, since /'(s) < on A' while ^'(s) > 0, we get that all s > so, 
(*,/(*))€ A'. 

On the other hand, since /'(s) < 0, for all s > si, /(s) < /(si) < 93(si). If x(ti) = s±, this 
in turn yields to : for all t > t\ , x"(t) = x'(t) — ip(x(t)) < /(si) — v?(si) = —5. Integrating, 
this implies that lim^oo x'(t) = —oo which contradicts x'(t) > 0. 

We have proved so far that for all t > 0, X(t) € A. This implies that x'(t) is increasing. 
In particular lim^oo x(t) = oo and S = oo. Also, x' being upper bounded by |M|oo> it must 
converge to a limit and lim^oo x"{t) = 0. Since 

lim <p(s) = IMloo, 
s— too 

we obtain the claimed statement. ■ 
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5.3 Comparison of second order differential equations 

For two functions ipi, ip2 in T^i we write ip\ < ip2 if for all t >0, ipi(t) < ip2{t)- 

Lemma 20. Let ip G Dq and x G T~L such that x = T v x and x'(0) > 0. Let ip G T> and y be a 
solution of y" — y' + ip(y) = with y(0) = 0, y'(0) > and. We define the exit times 

T = inf{t > : (y(t),y'(t)) £ 

T+ = inf{t > : y'{t) > |M|oo} and T_ = inf{t > : ip(y(t)) < y'{t)}. 

(i) IfT + < T, T + < oo and ip < ip then y'(0) > x'(0). 

(ii) IfT- < T, T_ < oo and ip > ip then x'(0) > y'(0). 

Proof: Let us start with the hypothesis of (ji|). The proof is by contradiction : we also assume 
that y'(0) < x'(0). We set Y{t) = (y(t),y'(t)) and G( yi ,y 2 ) = (y 2 ,V2 - ^(yi))- Define the 
trajectories $ = {X(t) eR^:t>0}, and for r > 0, *(r) = {Y(t) G R 2 + : < t < r}. By 
Lemma \W\ <3? is the graph of an increasing function / : M + —> R + with /(0) = x'(0) > : 

$ = {(*,/(*)) :seR + }. 

Similarly, if y G R+, G{y)\ = y2 > 0. Thus there exits a differentiable function g : [0, y(T)] — > 
such that 

*(T) = {( S)ff ( S )): S G[0 )2 /(r)]}, 

with 

f(s) 

Now, the assumption < y'(0) < x'(0) reads < g(0) < /(0). 

By assumption, there is a time S 1 > such that g(S) > ||</?||oo > f(S). Hence, by the 
intermediate value Theorem, there exists a first time < s\ < S such that the curves intersect 
: g(s\) = f(si) and g(s) < f(s) on [0, si). However, it follows from ([31)1 that for s G [0, si), 

9 '( S ) = 1-M <1 _^£| = m 

/(«) 

Hence, integrating over [0, si], we find <?(si) — g(0) = /(si) — <?(0) < /(si) — /(0). This 
contradicts g(0) < /(0). We have proved (0). The proof of ([u]) is identical and is omitted. ■ 



5.4 Proof of Theorem [7] 

We first linearize (|4U|) with y(s) = A(l — e _s ) in the neighborhood of A = 1/4. 
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Step one : Linearization from above. We have <//(0) = A, and from the concavity of ip, 

tp(s) < As. (45) 
We take A > 1/4 and consider the linearized differential equation 

y" - y' + Xy = 0. (46) 
The solutions of this differential equation are 

y(t) = asm(ut)e2 + 6cos(o;t)e2 , 

with 



jj ■ \ A - ^. 



We use this ODE to upper bound x' x (0). 
Lemma 21. For all A > 1/4, 

^ ( 0)<^ e -£(l + O(u, 2 )). 
Proof: Let a > and consider the function 

y{t) = asm(ujt)e2 . 

We have y(0) = 0, y'(0) = aw, 

y'(t) = ae2 (a; cos(a;t) + -sin(a;t)), 

and 

y"(i) = aea ^a;cos(a;t) + ^- — a; 2 ^ sin(u;£) 
Define 

T = — arctan [ = - - ) = - - 4 + 0(w 2 ). 

On the interval [0, T], y"(t) > and y"(T) = 0. Thus the function y'(t) is increasing on [0, T] 
and 

y '(T) = e~ 2 ae^(w + 0(a; 3 )). 
Hence, we may choose a such that y'(T) = A = I + w 2 with 

a = |^(l + 0( w 2 )). 
It remains to use (|45p with Lemma [20l ■ 
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Step two : linearization from below. For < k < 1, we define 

i = - + k 2 uj 2 < A, 
4 

and the function in T> 

ip(s) = min (Is, <p(s)) . 

In particular 

We shall thus consider the linear differential equation 

y" -y' + £y = 0, 

The solutions of ([35]) are 

y(t) = a sm(unt)e2 + bcos(uJKt)e2 . 

A careful choice of a, k will lead to the following lower bound. 
Lemma 22. For all A > 1/4, 

4(0) > — wV£(l + 0(w 2 )). 

7T 

Proof: For a > 0, we look at the solution 

y(t) = asin(a;Kt)e2. 

We have y(0) = 0, y'(0) = anu. 

y'(t) = ae? (uk cos (oont) + — sin(u)Kt)). 
We repeat the argument of Lemma I2T1 On the interval [0, T], y"(t) > and y"(T) = 0, 

r = arctan - — = 4 + 0(w 2 , 

and the O(-) is uniform over all k > 1/2. The function is increasing on [0, T] and 

y '(T) = ae~ 2 e^ujK(l + +0(oj 2 )). 
Now, we have is < y(s) for all s G [0, <r] with 

&r = A(l - e~ CT ). 

It yields to 

a = 2^1 - £) + of 1 - ~Y = 8(1 - k 2 ^)uj 2 + 0((1 - k 2 )w 4 )- 

30 



However from f|48j) . for t = T, since y"(T) = 0, we have 

y(T) 4 
We may choose a such that y(T) = a by setting 

a = at- "—(I + 0{u 2 )) = 2e"e-^^h(l + 0{u 2 )). 

4 K 

Now, in the domain < y < a the non-linear differential equation y" — y' + ^(y) obviously 
coincides with (j4~8l) . Thus, using (jlTj) and Lemma [20], it yields to 

x'(0) > y'(0) = 2eV^(l - k 2 )uj 2 {\ + 0(oj 2 )). 

Taking finally k = 1 — 2u/tt gives the statement. ■ 



Step three : End of proof. We now complete the proof of Theorem [3 We start with the 
left hand side inequality. By (I42p . we have 

x{t) > x'(0)t. 

It follows from (|36p that 

poo poo -1 

qW = l ^ e-<°>< rtt = _i_. 

It remains to use Lemma [22] and we obtain the left hand side of Theorem [7] 

We now turn the right hand side inequality. For X = {xx,Xq) 6 M 2 , define G(X) = 
(x2,X2)- From the definition of F in (J43]), we have, component- wise, for any X £ R 2 , 

F(X) < G(X). 

Note also that G is monotone : if component- wise X < Y then G(X) < G(Y). It follows 
easily that if X(0) = Y(Q), X' = F{X) and Y' = G{Y) then 

X(t) < Y(t). 

Looking at the solution of y" — y' = such that y(0) = and y'(0) = x'(0), we get that 

x(t) < x'(0)(e* - 1). 
We deduce from ([36]) that, for any T > 0, 

poo pT 

q(\) = / e-^e^dt > / e~ x ' ^ e~ l dt 
Jo Jo 

> [ (1 -x'(0)(e* - l))e-*dt 
•/ o 

> 1 - e _T - x'(0)T. 

Now, we notice that in order to prove Theorem [71 by Lemma [2~T| we may choose A close 
enough to 1/4 so that x'(0) < 1. We finally take T = — ln(x'(0)) and apply Lemma[2TJ This 
concludes the proof of Theorem [7j ■ 
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Appendix 



In this appendix, for the sake of completeness we include the proof of the following lemma on 
Galton- Watson trees. 

Lemma 23. Let T be a Galton-Watson tree with mean number of offsprings d > 1. Condi- 
tioned on T is infinite, T is a.s. lower d-ary. 

Proof: Let 1 < 5 < d and Z n = \V n \ be the number of offsprings of generation n. From 
Seneta-Heyde Theorem (see [20^ Chapter 5]), conditioned on T is infinite, a.s. 

lim — log Z n = d. 

n— >-oo n 

Let p > be the probability that T is infinite. It implies that for any e > 0, for all n large 
enough, we have ¥(Z n > 5") > p — e. 

Now, consider a new Galton-Watson tree T' starting from the root where each vertex 
produces independently m = [5 n \ offsprings with probability p — e and offspring otherwise. 
From what precedes, we may couple T and T' such that T" is a subtree of T* n . 

We are now going to prove that T' contains a large regular tree with probability at least 
p — 2e. To this end, we set 

q = q(e) = 1—p + s. 
Note that we may have chosen n = n(e) large enough so that 

q + (1 - q)e- m£2/2 <q + e. (49) 

We consider the following pruning algorithm on T'. At step 0, we start with all vertices of 
T' . At step 1, we remove all vertices which have less than (1 — q — 2e)m offsprings. We now 
iterate: at step k > 1, we remove all vertices which have less than (1 — q — 2e)m offsprings 
left by step k — 1. 

Denote by pk the probability that the root of T' is removed by step k. We have po = 0, 
pi = q and (pk)k>o is an non-decreasing sequence. We are going to check by recursion that 
for all k > 1, 

p k <q + e. (50) 

Indeed, let k > 1 and assume that Pk-i < q + £• Note that if the root is removed by step 
k, then either it has offspring or more than (q + 2e)m of its offsprings were removed by 
step k — 1. From the recursive structure of the Galton-Watson tree, the probability that an 
offspring was removed by step A; — 1 is pk~i and these events for each offspring are independent. 
It follows that 



Pk <q+(l-q)F \Y, X i ^ {q + 2e)m 



where (Xj)i<j< m are i.i.d. {0, l}-Bernoulli variables with mean pk-i- By recursion hypothesis, 
Pk-i < q + £• Hence, Hoeffding's inequality leads to 




Xi - (? + 2e ) m ) - p ( ~ ^ £m ) - e 



-me 1 12 
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From ([4"9~|) . we deduce that pk < q + e. This proves ([50]) . 

We have thus proven that with probability at least 1 — (q + e) = p — 2s, the root of T is 
never removed by the pruning algorithm. However, on the latter event, by construction T' 
contains a [(1 — q — 2e)mJ-ary tree rooted at (note that 1 — q — 2e = p — 3s). 

We may now conclude the proof. We apply the above argument to some 5' £ (5,d). This 
proves that for any < e < 1, there exists an integer n £ such that with probability at least 
(1 — e)p, T* ne contains a [~<5 ne ]-ary tree. Note that the latter event is contained in the event 
that T is infinite. It follows that the conditional probability that T* n£ contains a [<5 ne ]-ary 
tree, given T infinite, is at least 1 — s. 

We finally consider the sequence £k = 1/k 2 . From Borel-Cantelli lemma, conditioned on 
T infinite, a.s. there exists k such that T* n£ * contains a [~<5 n£ fc ] -ary tree. ■ 
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